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ABSTRACT 

We propose a new method for the study of the chiral properties of the 
ground state in QFT's based on the computation of the probabihty distri- 
bution function of the chiral condensate. It can be applied directly in the 
chiral limit and therefore no mass extrapolations are needed. Furthermore 
this approach allows to write up equations relating the chiral condensate with 
quantities computable by standard numerical methods, the functional form 
of these relations depending on the broken symmetry group. As a check, we 
report some results for the compact Schwinger model. 
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The study of the chiral properties of the ground state in Quantum Field 
Theories (QFT's) is one of the main points of attention in recent develop- 
ments of Lattice Gauge Theories (LGT's). In QCD, the gauge theory for the 
strong interacting sector of the Standard Model, the low energy physics could 
not be understood without a detailed analysis of the chiral properties of the 
vacuum. Furthermore there are strong evidences from numerical simulations 
of QCD at finite temperature, suggesting that a chiral transition appears at 
some critical temperature Tc separating the quark-gluon plasma phase from 
the hadronic phase. 

Concerning the abelian model, recent numerical studies in the noncom- 
pact formulation show the existence of a strongly coupled QED and open 
the possibility to get a non-trivial continuum limit for a non asymptotically 
free gauge theory in four dimensions. The strongly coupled phase of this 
model is characterized by a vacuum state not invariant under chiral transfor- 
mations and therefore a mechanism for dynamical mass generation appears. 

Last, the analysis of the chiral properties of the vacuum in 2+1 QFT's and 
of their dependence on the flavor number has also become of great interest be- 
cause of its possible relevance in the high Tc superconductivity phenomenon. 

In this paper we propose a new method for the study of the chiral prop- 
erties of the ground state in QFT's which is based on the computation of 
the probability distribution function (p.d.f.) of the chiral condensate in the 
chiral limit. This is an standard procedure when analysing spontaneous sym- 
metry breaking in spin systems or in QFT's with bosonic degrees of freedom, 
since this kind of degrees of freedom can be simulated directly in a computer. 
However, in in the numerical simulations of a LGT with dynamical fermions, 
the Grassmann fields must be integrated analytically for obvious reasons. 
Then, even if the ground state is degenerate, in the analytical procedure we 
integrate over all possible vacuum states, the chiral order parameter being 
always zero for massless fermions. We will show here how despite of the fact 
that Grassmann variables cannot be simulated in a computer, an analysis 
of spontaneous symmetry breaking without a symmetry breaking external 
field, which is standard in te case of spin systems, can also be done in QFT's 
with fermion degrees of freedom. Previous attemps to analyse the vacuum 
structure in fermion-gauge systems |^ were based on several approximations. 
We want to remark that our analysis is completely free from approximations. 

The main advantage of this method, when compared with standard simu- 
lations, is that we can work directly in the chiral limit and therefore no mass 
extrapolations are needed. Furthermore this approach allows to write up 
equations relating the chiral condensate with quantities computable by stan- 
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dard numerical methods, the functional form of these relations depending on 
the nature of the broken symmetry group. 

Another interesting field of application of p.d.f. of order parameters is 
the quantitative description of interacting social groups in stochastic models 
for the formation of Public Opinion. Previous work in this field analyses 
the time dependent p.d.f. of public opinion for the simplest case of two kind 
of opinions (+, — ), modeling it with an Ising ferromagnet, the dynamics of 
which takes into account the influence on an individual of its neighbors. 
Generalization of this work to the case of a continuous opinion (continuous 
symmetry) could be another possible application of the formalism developed 
here. The application of fermionic models to Public Opinion formation is 
interesting not only for the continuous character of the symmetry but also 
because of the long range interactions induced by the fermion dynamics, 
which would give a more precise description of the high degree of long distance 
communication in present world. 

1. Theoretical Grounds 

Our starting point is a QFT describing a gauge field coupled to a fermion 
matter field and regularized by means of a space-time lattice. Next suppose 
the ground or equilibrium state of this model is degenerate, and be a the 
index which characterize all possible vacuum states. The ensemble of all 
equilibrium states is a partition of the Gibbs state and the probability Wa 
to get the vacuum state a when choosing randomly an equilibrium state 
depends on the total free energy of the a state . 

Since we are interested in the analysis of the chiral properties of the 
ground state, let us choose as order parameter the chiral condensate ipip and 
characterize each vacuum state a by the expectation value Cq of the order 
parameter in the a state, i.e. 



where N is the total number of lattice sites and the sum is over all lattice 
points. The p.d.f. P{c) of the chiral order parameter c will be given by 
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The function P{c) tells us what is the probability that choosing randomly 
a vacuum state, we get the value c for the chiral order parameter. If the 
vacuum state is invariant under chiral transformations, i.e., if it is unique 
as concerning the chiral symmetry, P{c) will be a single S function S{c). 
Otherwise P{c) will be a more complex function, sum of 6 functions in the 
case of a discrete symmetry group or a continuous function in the other cases. 

The next step now is to relate -P(c) with magnitudes which can be com- 
puted by numerical simulations. To this end, let us write 

P(c) = Zim^^oo < -c) > (3) 

where the expectation value in (|]) is computed in the Gibbs state and the 
integration measure is that associated to the partition function 

Z = J dipd'iljdUe-^°^^'^+'^^'^. (4) 

Sg in is the pure gauge action and A the fermionic matrix. 

To check that equation (|^) gives us indeed the p.d.f. of the chiral con- 
densate, it is enough to verify that all the moments of the p.d.f. defined in 
(H) agree with those of (|^). The verification for the first moment is triv- 
ial whereas for the higher moments it is enough to take into account that 
in the thermodynamical limit intensive quantities do not fluctuate in any 
equilibrium state a, i.e. 

Xl,X2,...Xp 

= ^ < 1piXi)ij{Xi) >„< ^(X2)^(X2) >a < ^l'{Xp)'lp{Xp) . (5) 

Xl,X2,...Xp 

Expression is not suitable for numerical computation. However we 
can define its Fourier transformed P{q) 

P{q) = J e''^'P{c)dc (6) 
which, as will be shown, can be numerically computed. 
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Since we are interested in the analysis of the chiral symmetry on the 
lattice, we will use the staggered fermions regularization. The fermion matrix 
A can be written as 



A = m + iA (7) 

where m is the fermion mass and A a hermitian matrix which depends on 
the gauge field configuration. The eigenvalues of A are real and symmetric. 
Taking into account all these properties of A, the following expression for 
P{q) can be derived from (|]) 



and in the chiral limit 



m =< n(i - > (9) 

where the product in (H) and @ runs over all positive eigenvalues Xj and 
the mean values are computed with the probability distribution function of 
the effective gauge theory obtained after integrating out the fermion fields. 
In order to get meaningful results for (P), special care must be taken with 
the boundary conditions for the fermion field to avoid zero modes. 

The function P{q) can be computed numerically and then, by inverse 
Fourier transform we get -P(c). We can go deeper in the investigation of 
the form of the p.d.f. P{c). In the lattice regularized action of a gauge the- 
ory with Kogut-Susskind fermions, only a U{1) subgroup of the continuous 
chiral symmetry is preserved. Then if this continuous residual symmetry is 
spontaneously broken, we will get a continuum of equilibrium states charac- 
terized by an angle a with a taking values between — vr and vr. The v.e.v. of 
the chiral condensate at each vacuum will be given by Cocos(2a), cq being 
the value corresponding to the a — vacuum selected when switching-on an 
external "magnetic" field. Therefore, the function P{c) can be computed as 

P(c) = ^ r daS{c - Co cos(2a)) (10) 
2tt J-tt 
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which gives for P(c) the value l/(7r(cQ — c^Y^'^) for — cq < c < cq, P(c) = 
otherwise (see Fig. 1). Its Fourier transformed 



P{q) = — r dOe'"'"''''^ (11) 
27r J-TT 

is the well known zeroth order Bessel function of the first kind Jo{qco). 

Several relations between the chiral condensate and the eigenvalues of the 
fermionic matrix can be derived from (|l5). For example, by doing the second 
derivative of the function P{q) we get the second moment of the distribution 
P(c) and then the following relation holds 



co' = <^^>' = <4E^>- (12) 

3 i 

In the symmetric phase, the right hand side of (0) is the longitudinal sus- 
ceptibility normalized by the lattice volume A^. Furthermore the function 
P{q) is proportional to the partition function Z evaluated in the imaginary 
axis of the complex mass plane and therefore a relation between the zeroes of 
the partition function in the complex mass plane and the value of the chiral 
condensate follows from this relation. 

Using ( pTD the following relations between the chiral condensate cq, the 
zeroes qi of P{q) and the zeroes of the zeroth order Bessel function of first 
kind can be derived. 



co(^) = -. (13) 

qi 

In the thermodynamical limit the value of Cq will be independent of the 
order i of the zero in (0). At finite volume Cq depends on i. However we 
have observed that a plateau appears in the plot of cq as a function of i 
when chiral symmetry is spontaneously broken. The extent of this plateau 
increases with the lattice size becoming eventually infinite for infinite lattice 
volumes. These things are very well illustrated in Fig. 2 where we plot some 
preliminary results on the lattice chiral condensate co{i) against i for the 
one-fiavor compact Schwinger model, obtained from MFA simulations [^]. 
The continuous line in this figure stands for the continuum analytical result 
(in units of charge) times Of course the Schwinger model is not very 

interesting from a physical point of view but, as well known, it is a very 
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good laboratory to check new proposals for two reasons: exact analytical 
results are available and it shares many interesting properties with other 
more relevant physical models. This is the reason why we believe the results 
reported in Fig. 2 are very encouraging. 



2. Discrete versus continuous p.d.f. 

The analysis of p.d.f. of order parameters in Statistical Mechanics is not 
new. It has been developed for ferromagnetic Ising-like systems and also 
for spin-glasses |^ , the last with a much more complicated vacuum structure 
than magnets. However this is the first time to our knowledge this kind of 
analysis is apphed, without any approximation, to QFT's with dynamical 
fermion fields. 

There are several important differences between p.d.f .'s in spin systems 
and fermionic systems and we would like to point out them here. First and 
from a practical point of view, the p.d.f. of the density of magnetization in 
a spin system can be directly measured by computer simulations whereas 
in a fermionic system this is not possible since Grassmann variables cannot 
be simulated in a computer. This is the reason why in our case we work 
with the Fourier transform P{q) rather than with P{c) (^). Furthermore 
and as a consequence of the anticommuting character of fermionic fields, 
all the moments of the p.d.f. P{c) of order larger than the lattice volume 
N vanish since they are correlation functions of anticommuting variables 
with some repeated index. In fact the function P{q) given by equation (P) 
is a polynomial of degree the lattice volume and therefore its inverse 
Fourier transform does not exists as a regular function if P{q) is defined in 
an unbounded momentum space. 

This peculiarity of the p.d.f. for fermionic systems complicates the finite 
size scaling analysis of it since in order to do such a kind of analysis, we need 
to put an arbitrary cut-off in momentum space which could induce extra 
finite-volume effects. However we can by-pass this difficulty by applying the 
finite size scaling hypothesis, instead to the p.d.f. p{c) as in Binder's work on 
Ising systems [^, to its Fourier transform P{q). Then the finite size scaling 
hypothesis for P{q) reads 



Pdq) 



(14) 



where L is the linear dimension of the lattice and / is a universal scaling 
function. Equation should work for L large and near the critical point. 
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The next step now should be to estimate the exphcit form of the universal 
function / in both the symmetric and broken phases. Since this kind of 
analysis is too long to be reported here, we will stress only some important 
differences which appear when comparing with Binder's analysis of Ising 
like systems. The physical origin of the different behavior of the universal 
function / for Ising and fermion-gauge systems is the existence of a continuum 
of equilibrium states in the broken phase of fermion-gauge systems in contrast 
to the two ground states of the low temperature phase of the Ising model. 

The symmetric phase of both models is characterized by a non degener- 
ate vacuum and Binder's analysis 0, which is based on the assumption that 
in this case -Pl(c) (equivalently Pl{(i)) is well approximated by a gaussian 
function centered at the origin, should work also for fermionic systems, at 
least for not too large values of q. The situation however changes drastically 
in the broken phase where Binder's analysis, based on the two gaussian ap- 
proximation of the p.d.f. at finite lattice size, fails completely to describe the 
spontaneous breaking of a continuous symmetry. In fact the infinite volume 
limit of Pl{<i) (see eq. ([Tl|)) is a Bessel function in contrast with the cosine 
function which appears in the Ising case. Taking into account this result 
and writing the universal function / in the broken phase as a function of the 
chiral condensate 



f{qL=^,j)=nqL=^,\)- (15) 
the following expression for the universal function F{z^ z') holds 

F{z,z') = Jo{^) + , (16) 

z 

where Jq is the Bessel function and the dots in ( [T6| ) stand for size dependent 
terms which vanish in the infinite volume limit. There are many ways to 
parameterize the first corrections in (|TB|), the simplest one being a term 
proportional to as follows from the fact that the universal function F and 
the Bessel function Jq differ in an analytical function of g^, which vanish at 
the origin in momentum space. A detailed analysis in this direction for some 
specific models will be reported in a separate paper. There are however some 
interesting features concerning finite volume effects which can be derived from 
the results reported in Fig. 2, and which we would like to point out. 
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Our analysis of the Schwinger model shows that, as expected, finite vol- 
ume effects are more or less relevant depending on the quantity computed. 
By looking at Fig. 2 we can understand why finite size effects are relatively 
large when computing moments of -P(c), like equation (p!2D. In fact in this 
case all the zeroes of P{q), or equivalently all the eigenvalues of the fermionic 
matrix, give contribution to the moments of -P(c) and our results show clearly 
that higher order zeroes of P{q) suffer from stronger finite size effects. What 
we find certainly surprising in the results of Fig. 2 is the fact that an ac- 
curate estimation of the chiral condensate in the Schwinger model can be 
obtained from rather small lattices, by analyzing the first zeroes of P{q)- 
They show an asymptotic scaling behavior which follows from the fact that 
the first eigenvalues of the dominant gauge field configurations scale with the 
inverse lattice volume, the window for this scaling increasing with the lattice 
size. These results should stimulate people working in this field to apply this 
formalism to more interesting physical systems, like QCD. 

We thank CICYT (Spain) for partial financial support to this work. 
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Figure captions 

Figure 1. Standard form of P{c) (cq. (1.11)) in the broken phase. 

Figure 2. Chiral order parameter in the compact Schwinger model in 
32^ and 64^ lattices at /5 = 6.344. The solid line corresponds to its continuum 
analytical value. 
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